QCD sum rules for quark-gluon three-body components in the B meson 
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We discuss the QCD sum rule calculation of the heavy-quark effective theory parameters, \e 
and Ah, which correspond to matrix elements representing quark-gluon three-body components in 
the B-meson wavefunction. We derive the sum rules for \e,h calculating the new higher-order 
QCD corrections, i.e., the order as radiative corrections to the Wilson coefficients associated with 
the dimension-5 quark-gluon mixed condensates, and the power corrections due to the dimension-6 
vacuum condensates. We find that the new radiative corrections significantly improve the stability 
of the corresponding Borel sum rules and lead to the reduction of the values of \e,h- We also 
discuss the renormalization-group improvement for the sum rules and present update on the values 
of \e,h- 
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I. INTRODUCTION 

The B mesons play distinguished roles in exploring CP 
violation and the flavor sector of the Standard Model. In 
particular, the measurements of the i?-meson decays can 
provide precise information on the relevant quark cou- 
plings [1]. Since the properties of those decays are also 
influenced by the complicated strong-interaction effects 
responsible for forming the B-meson, as well as the final- 
state hadrons, a better theoretical control of the nonper- 
turbative effects inside the B mesons is now becoming 
very important. This is also interesting in its own right 
as understanding the properties of the simplest meson 
including a heavy quark. 

In the heavy-quark limit based on Aqcd/tob ^ 1, 
AqcD/mf, ^ 1 with ms and being the masses of 
the i3-meson and fo-quark, the matrix elements using a 
B-meson state obey heavy-quark symmetry, and are con- 
veniently described by the heavy-quark effective theory 
(HQET) (21. In this framework, fundamental properties 
of the B mesons are represented by the HQET parame- 
ters that are defined as matrix elements of the relevant 
local operators, like the decay constant F [2]: 



\qiplbhv\B{v)) ^iF{fi)vp 



(1) 



Here, \B{v)) is the _B-meson state with the 4- velocity v in 
the HQET, q is the light-antiquark field, hv is the effec- 
tive heavy-quark field, and the heavy-light local operator 
in the LHS is renormalized at the scale fj,. The decay con- 
stant of Eq.([T|) represents the quantitative content of the 
quark-antiquark valence component inside the B meson 
in the heavy-quark limit, so that F(fj,) determines the 
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normalization of the valence Fock components in the B- 
meson wavefunction, as well as of the amplitude for the 
exclusive B-meson decays. We note that F{^) is related 
to the physical decay constant fs as 



1 



0{l/mb), 



CpCts 

47r 



31n^ 



(2) 



with the corresponding short-distance coefficient shown 
in the parentheses to the one-loop accuracy, as well as 
with the 0(l/mfc) correction terms in the heavy-quark 
expansion. The value of fs is now obtained rather pre- 
cisely from lattice QCD calculations as Q (see also Q) 



/s =0.195 ±0.013 GeV, 



(3) 



which is consistent Q with the results of measurement of 
the branching fraction for B ^ tv decays in the Belle @ 
and Babar [7| experiments. 

We can also define the analogues of Eq. ([Ij , which are 
associated with the higher Fock components inside the B 
meson. For the non-minimal parton configurations with 
additional gluons, the corresponding HQET parameters 
were introduced by Grozin and Neubert [Sj as 

mac ■ gEj5K\B{v)) = F{p)Xl{^^), (4) 
{0\qcT ■ gH-izK\B{v)) = lF{^i)\l{^i), (5) 

in terms of the matrix elements in the B-meson rest frame 
with V — (1,0). Here, the three-body quark-gluon opera- 
tors are associated with the chromoelectric and chromo- 
magnetic fields, = and B' = (-1/2)£'J'=G^'=, with 
G^. = G^,r- and G^, = d^A- ~ d^A- + gf^'^^A^Al 
being the gluon field strength tensor. The values of \\. ^ 
were also estimated in Q using QCD sum rules, as 
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=0.11 ±0.06 GeV^ 
A|^(/x) = 0.18 ±0.07 GeV^ 



(6) 
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at /Lt = 1 GeV. Besides this rather rough estimate, there 
exists no other estimate at present. In this paper, we 
present an extension of Grozin-Neubert's QCD sum rule 
calculation of A|; ^, taking into account the higher-order 
perturbative and nonperturbative effects in QCD. The 
main new ingredient in the present case is that we cal- 
culate the relevant QCD radiative corrections, so that 
we derive our sum rules to the order accuracy. We 
find that, only after including those new contributions, 
the perturbative as well as nonperturbative corrections to 
the sum rules for X"^ ^ become under control. We note 
that, also for the QCD sum-rule calculations of the decay 
constant ([ij, the order Ug radiative corrections produce 
the large and essential contributions to yield the values 
consistent with Eq.® (3,i,[i3l- 

One might anticipate that the higher Fock components 
in the B meson would give rise to the "higher-twist" 
power corrections to the hard exclusive amplitudes, sim- 
ilarly as the roles played by the higher Fock components 
in the light mesons tt, p, etc. jTl| . and thus the impact 
of a more accurate determination of \% fj than Eq.([6]) 
would be marginal. Actually, however, it has been re- 
vealed that the behaviors of the contributions induced 
by the higher Fock components are quite different be- 
tween the B-meson case and the light-meson case: the 
presence of a heavy quark inside the _B-meson causes 
the nonperturbative quark-gluon interactions which in- 
duce the mixin g of the effects corresponding to the dif- 
ferent twist [sl. Il2l. ITs! ]. In particular, recently, it has 
been demonstrated that the B-meson "light-cone distri- 
bution amplitudes" to describe the valence Fock conipo- 
nents participating in the hard exclusive processes [1, u3r - 
[TgI are contaminated by the multiparticle Fock states, so 
that the contributions represented by the novel HQET 
parameters A|, ^ of Eqs. ([U, ([5]) could strongly affect [13] 
the amplitudes for the exclusive i?-meson decays at the 
leading power. Indeed, the normalization of the so-called 
hard spectator interaction amplitude [1[ could be mod- 
ified by a factor two or more, when varyin g th e values 
of Xi, jj in the uncertainty range ofEq.® 17]. There- 
fore, an improved estimate of \\ is desirable to have 
a better control of the hadronic uncertainty associated 
with the B meson, which is a major source of theoretical 
uncertainty in the calculations of the decay rates [l|, [3] ■ 

The paper is organized as follows. Section |ll] is mainly 
introductory; we set up a systematic formalism for the 
QCD sum rule calculation of \\ ^ in the HQET, as well 
as the decay constant F, and apply it to reproduce the 
previous results for the sum rules of those HQET param- 
eters. We explain that the previous sum-rule estimate 
of X\ ^ needs update including higher-order effects. In 
Sec. Ilill we derive the new power corrections to the sum 
rules for A|^^, due to the nonperturbative QCD conden- 
sates. Section|lV]is devoted to the calculation of the new 
order as corrections to the sum rules for Xj^ ^ . Taking 
into account all these new contributions, we present the 
final form of our sum rules for Xj^ ^ in Sec. |Vl We ex- 
plain the renormalization-group improvement of our sum 



rule formulas and perform their detailed numerical anal- 
ysis to obtain a new estimate of A|; ^. We find that the 
new values of A|; are significantly modified from those 
of Eq.®. Section IVll is reserved for conclusions. 

II. QCD SUM RULES IN THE HQET 

The aim of the present paper is to update quantitative 
estimates of the S-meson matrix elements A_b and A^. 
To this end we apply the QCD sum rule approach to the 
suitable correlation functions in the HQET. In particu- 
lar, in this section, we set up the framework convenient 
for treating the perturbative as well as nonperturbative 
corrections to the relevant correlation functions and also 
use it to demonstrate that the calculation of the leading 
effects reproduces the sum rules obtained previously by 
Grozin and Neubert Q . The complete treatment includ- 
ing the new perturbative and nonperturbative corrections 
is presented in the succeeding sections. 

We consider the following correlation functions in the 
HQET (the dependence on the renormalization scale fi is 
suppressed for simplicity): 

i J d^xe-'"^-^{0\T[qiO)'r,h.,{0)K{x)T2qix)] |0) 

= -iTr[riP+r2]H;^H, (7) 

i J d^xe-^'^^-'^ {0\T [q{0)rigG^,{Q)h,{Q)h,{x)r2q{x)] |0) 

= -iTr[a^,riF+r2]H3ff(w) 

-^Tr [(wp7^ - iv^j^)TiP+T2] H3s(a;), (8) 

where P+ = (1 -I- 1^)/2 is the projector on the upper com- 
ponents of the heavy-quark spinor, Fi is an arbitrary 
gamma matrix, and we choose F2 — 75 to construct 
the sum rules for pseudoscalar B meson. (The case for 
the vector meson B* can also be treated by choosing 
^2 = ~ and yields exactly the same results as in 
the pseudoscalar i?-meson case due to heavy-quark spin 
symmetry in the HQET.) Eq. ([7]) is the familiar correla- 
tor between two heavy-light currents and the correlation 
function Hpi'^j) provides the sum rules to evaluate the 
decay constant F. On the other hand, the correlator 
between the two-body current and the three-body cur- 
rent involving the gluon field strength tensor defines the 
correlation functions, \\^h{lo) and H35(cli), correspond- 
ing to the two independent Lorentz structures; as we will 
show shortly, W^h{'^) and H35(a;) allow us to derive the 
sum rules to evaluate A|^ and the "splitting" A|^ — 
respectively (see Eqs. (jlj, ([5])). 

In the general procedure of QCD sum rules, we evalu- 
ate the above correlation functions by the operator prod- 
uct expansion (OPE) in the unphysical region —u ^ 
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Aqcd on one hand and express those correlation func- 
tions in terms of the properties (masses and matrix ele- 
ments) associated with the physical states participating 
in the spectra at w > on the other hand; we relate 
these two descriptions exploiting the analyticity proper- 
ties of the correlation functions, which are embodied by 
the corresponding dispersion relations. The dispersion 
relation satisfied by Ilp{uj) of Eq. (O is well-known, and 
the dispersion relations of similar type are obeyed also by 
n3jy(w) and n3s(w) of ©; namely, for X = F, 3H, 35, 

, Imnx(w') 



w' — w — iO' 



(9) 



up to the appropriate subtraction terms that are polyno- 
mial in uj. Indeed, these relations can be demonstrated 
by inserting a complete set of states between the two cur- 
rents in the corresponding correlators: the LHS of Eq. ([7]) 
yields (in the S-meson rest frame) 

^ '~\qmihAO)\B{v)) 



2{A-uj-iO) 

x{B{v)\h,{0)T2qmO) 



(10) 



and the LHS of Eq. ([8]) gives the similar result with Fi re- 
placed by ri(7G'^i^(0). Here, a pole arises at A = ms— m;,, 
the usual effective mass of the B meson, accompanying 
the matrix elements that are parameterized by the cor- 
responding HQET parameters (renormalized at the scale 
as 



B{v)) ^--F{fi)TT[T,P+j5], (11) 



(0|g(0)ri/i, 

= -^F(/.){Al,(M)Tr[riF+75fT^.] 

and the ellipses stand for the similar pole contributions 
of higher resonances and the continuum contributions. 
Combining the results with Eqs. ([7]), ([S]), one finds, 



n3_ff(w) = 
n3s('^) = 



1 



2 A-uj-iO 
F2(^).,_ 1 
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A 



iO 



A- uj - iO 



•(A3) 



which hold in any frame as well as in the rest frame, 
and show that the contribution of the B meson to the 
spectral functions in Eq. ([9]) is completely expressed by 
the relevant HQET parameters. Calculating the LHS of 
Eqs. (O, ([S]) based on the OPE and matching the results 
with the formulas in Eq. ljlSp . we obtain the sum rules 
associated with those HQET parameters. 

We follow the standard procedure to construct the 
corresponding QCD sum rules: we apply the Borel- 
transformation operator, defined by 



Bm — 



lim 

AI——uj/n fixed 



Tin) 



d_ 

duj 



(14) 



to the relevant correlation functions obeying the disper- 
sion relation This transformation introduces 
the Borel parameter M instead of the external energy ui 
as 



1 f°° 1 

' dw'e-"/^^-ImHx(w'), (15) 

TT 



and eliminates the subtraction terms. Eg. (1141) implies 
that the power-correction terms associated with the 
higher dimensional operators in the OPE are factorially 
(~ 1 /n!) suppressed, improving the convergence of the se- 
ries and, simultaneously, Eq.([T5]) indicates that the con- 
tributions of higher resonances and continuum are expo- 
nentially suppressed compared with that of the lowest- 
lying state. In particular, owing to the latter property 
of the Borel transform to minimize the dependence on 
the contributions of the excited states, we approximate, 
as usual, those excited-state contributions to the spec- 
tral function in Eq. ljlSp by the continuum contribution 
which is based on the OPE result and starts from the 
"continuum threshold" Wth! namely, we use 



1 / N 

— Imni?(w) 

TT 



A) 



(16) 



with the correlation function H°^^(w) calculated in the 
OPE, and the similar form for (l/7r)ImHx(w) {X — 
3H,3S) with the corresponding OPE result, n^^^{uj). 
Then, we obtain the sum rules. 



F^^i)e-^/^' = 2 / dc.e-'^/^^-ImHOPi=(a 

Jo 



(17) 



Jo 



F\^.) [xl{^.)-xU^^)] e-'^"^ 

Jo 



(19) 



and, taking the ratios of Eqs. ([T8|) and ([T9|) with Eq. lfTT)) 
to cancel the factor F'^{fi), one can evaluate A|^(/i) and 
A|^(/i) — A|;(/i), respectively, based on the sum rules. 

Now the remaining task is to calculate H^^^(a;), to be 
substituted into the RHS of Eqs. (IT7l)-(I19l). Carrying out 
the OPE of the corresponding correlation functions for 
the region — w ^ Aqcd, the results generically take the 
form. 



H 



OPE/ 



X (c) = Cf (c.) + {L0){qq) + {u){G') 



C^{u){qgG-aq) 



(20) 



with X = F,3H,3S, where (qq) = {0\qq\0), (G^) = 
(0|(G^j2|o)^ and {qgG ■ aq) = maG^'" a ^,,q\Q) are the 
quark condensate, the gluon condensate, and the quark- 
gluon-mixed condensate, respectively, as the vacuum ex- 
pectation values of the dimension-3, -4, and -5 local oper- 
ators, and are associated with the corresponding Wilson 
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coefficients C^{lS) with k = q,G, and a; an increase in 
dimension of the operators imphes extra powers of 1/w 
for the corresponding Wilson coefficients, and the eUipses 
in Eq. ipni denote the terms with the operators of dimen- 
sion d > 6. Cf{uj), associated with the unit operator, 
coincides with the purely perturbative contribution to 
Hxi^j). The condensates as well as the coefficient func- 
tions in general depend on the renormalization scale /i. 

For the correlation function (O with the two-body cur- 
rents, the OPE can be derived in a standard way and the 
Wilson coefficients appearing in Eq. ([^0)) with X — F are 
obtained d, [13, m, is HI by evaluating the familiar 
Feynman diagrams, which involve the heavy-quark prop- 
agator in a background gluon field A^, 

I I 

xVexpfig f dsv ■ A{sv)] , (21) 
\ <y v-x / 

for the case of the D dimensions, where x^^ ^ — {v ■ 
x)v^ ^ and the factor in the second line, i.e., the straight 
Wilson line along the velocity v with the path-ordering 
operator allows us to organize the interactions with 
the gluon field exactly in the HQET. Thus, it is con- 
venient to use the Fock-Schwinger gauge, x'^ A^J,(x) = 0, 
for the background gluon field, so that the heavy quark 
does not interact with the nonperturbative gluons in the 
calculation for power corrections to the correlation func- 
tion ([7]) (see Fig. [^Ja) below). In this case, it is also 
well-known that we have very useful relations [2l|: for 
the classical background gluon field, 



duux^ Gini(ux\ 



(22) 



and, for the light-quark propagator, 
I I 27r^ (-x2 + iO) 2 



iT 



' D 

K 1 



1) 



327r-¥ (- 



-{^,a'-}5G^,(0) + ... ,(23) 



with the ellipses denoting the terms associated with op- 
erators of dimension d > 3; it is worth noting that the 
term associated with the dimension-4 operator is ab- 
sent from the ellipses. We do not give the details of the 
calculation of the Wilson coefficients here but, for later 
convenience, sketch the relevant steps: we decompose the 
quark as well as gluon fields into the "quantum" and 
"classical" parts; the quantum parts are contracted to 
yield the propagators in the classical background fields, 
like Eqs. (|2T|) . ([23)) . while the classical parts satisfy the 
classical equations of motion, = 0, ?; • = 0, and 
-C'GJ^i/ = g^qi q'lnT°'l' with the summation over all 
quark fiavors. For the matching at the leading accuracy 



in as, the correlator in the LHS of Eq.(l7]) is evaluated as 

T [g(0)ri/i„(0) K{x)T2q{x)\ = q{Q)TiK{Q)K{x)T2q{x) 

I ^ I 

+g(0)ri/i„(0)7^,(a:)r29(a;), (24) 
I I 



and, substituting Eqs. ([2T |) - (|23p into the first term in the 
RHS, we immediately obtain Cf in Eg. ipH)) as the purely 
perturbative contribution and also find that Cq vanishes 
up to the corrections of 0{a1), as a direct consequence of 
the absence of the operator G^ in Ea. (l23t as noted above. 
On the other hand, the vacuum expectation value of the 
second term in the RHS of Eg. ([Ml) yields, 

e{-v ■ x)5(3) (3.^) (0|g(0)riP+r29(a;)|0), 

= e{-v ■ x)5(3) {xj_) [(o|g(o)riP+r2g(o)|o) 
+x'^(o|g(o)riP+r2i?^g(o)|o) 

.M^'^/n|5(o)riP+r2i?^i?.(z(o)|o) + • • • ] 

= 0(-i; • x)5(3) (a;^)Tr[riP+r2] 



1 

-- X''X' 

2 



16 



X {qgG ■ aq) 



(25) 



where we used the equations of motion Ipq = in the 
last equality, and this allows us to obtain G^ as well 
as Cj. As the result, the relevant Wilson coefficients 
read [i,[i3,[2| 



riF( ^ 2, 

C, {.) ^ In— 

C^{u) = 0, 



16w3' 



(26) 



up to the terms polynomial in uj. We calculate the dis- 
continuities of Eq. ([20|) with these coefficients, across the 
cut along the line > in the complex uj plane, and sub- 
stitute the results into the RHS of Eq.dlTl). This yields 
the sum rule. 



1 



M 



16M2 



(qgG-aq), 



where the function, 

W'-"'\x) = 1 



EX' 
T 



k=0 



k\ 



(27) 



(28) 



arises from the integral over the duality region, < w < 
wth, in Eg. pT)) . This sum rule can be used for a leading 
estimate of the decay constant F{iJ,). 

The sum rules for a leading estimate of Xf^ ^ can be 
derived similarly. The corresponding calculation was 
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FIG. 1: Nonvanishing diagrams for the OPE of the correlation 
function ((8| in the Fock-Schwinger gauge. The double line 
denotes the propagator of the effective heavy quark, and the 
circle and box represent the two- and three-body currents, 
respectively. The four diagrams (a)-(d) generate, respectively, 
the first four terms of (HOj with X = 3H, 3S. 



(a) 



(b) 



FIG. 2: The vanishing subdiagrams. (a) is associated with 
classical background gluon field of Ea. (|22p . (b) and (c) have 
a loop involving quantum gluon field emanating from the field 
strength tensor in the three-body current. The external lines 
with a bar at their end are amputated. 



performed by Grozin and Neubert [8| with two partic- 
ular choices for the gamma matrix Fi of ([8]), which 
make the corresponding three-body currents coincide 
with the chromoelectric and chromomagnetic operators 
in the LHS of Eqs. ^ and ([S]). For later convenience and 
for a cross-check of Grozin-Neubert's result, we here per- 
form the corresponding calculation for arbitrary gamma 
matrix Fi, and summarize the procedures and the re- 
sults. The correlator between the two-body and three- 
body currents in the LHS of Eq. dS]) is evaluated as 



T [g(0)Fi5G^,(0)/i„(0) /i,(x)F2g(a;)] 
= q{0)rigG^40)h40)h4x)r2q{x) 

I I 1 I 



^g(0)Fi.gG^,(0) KiO)Kix) T2qix) 



(29) 



By contrast to the above case leading to the results ([26|) , 
the extra gluons emanating from the three-body current 
participate in the present case. Those extra gluons can 
interact with the light-quark and such contributions re- 
quire the participation of the additional quark-gluon cou- 
pling in perturbation theory, so as to form the propaga- 



tor. 



I I 



F (f ) S^" 



+- 



F (f - 1) r 



27r^ (-z2 +iO)~ 

be 

(30) 



Stt- (-z2 -hiO) 
^gGlp{Q)zPg^x + '^gGl^z^ 



for the case of the D dimensions and using the Feynman 
gauge for the quantum part of the gluon field, with the 
ellipses denoting the terms associated with operators of 
dimension d > 3. The ellipses in Eq.(|29|) stand for the 
terms of this type at the leading accuracy in as, which 
are induced by the first term in the RHS of Eq. (l30l) , and 
the corresponding nonvanishing contributions are repre- 
sented by the Feynman diagrams (a) and (b) in Fig. [T] 
note that, by explicit calculation, the subdiagrams (b), 
(c) in Fig. [2] vanish, reflecting that a physical gluon repre- 
sented by the field strength tensor does not interact with 
the heavy quark, nor is absorbed into a single quark on 
the mass shell. (The contributions induced by the second 
term of Ea. (pO|) will be discussed in Sec. lIVn Decompos- 
ing those contributions from Figs. [T] (a) and (b) into inde- 
pendent Lorentz structures, as in the RHS of Eq.®, we 
obtain the Wilson coefficients and , respectively. 



in the OPE ^ with X = 3iJ, iS. 

Similarly, it is straightforward to see that the vacuum 
expectation value of the first term in the RHS of Eq. ip^ 
yields Gq '' , by combining the field strength tensor of 
the first term of Eq. ip^ with that from the second term 
of the quark propagator (j23p . corresponding to the dia- 
gram (c) in Fig.[TJ On the other hand, the contribution of 
the second term in the RHS of Eq. ip^ can be evaluated 
similarly as Eq.(|25|). and we obtain. 



e{-v ■ x)5(3) (a;^) (0|g(0)gG^,(0)FiP+r2(z(a:)|0) 



X I — Tr [FiP+Fsa^.] {qgG ■ aq) 



(31) 



with the term represented by the diagram (d) in Fig. [T] 
and the ellipses denoting the contributions associated 
with the operators of dimension d > 6, so that we can 
calculate G^^'^^ using the former contribution. Collect- 
ing the results from the diagrams (a)-(d) in Fig[l] we 
obtain the corresponding Wilson coefficients as 



and 



Gnu) = 

G%"{u:) = 



In- 



187r2 47r 

Cp — win , 

27r ij, 

247r " ^ ' 
1 

2Av' 



J^NcGp—uj^ln- 
187r2 ^ 4tt 



(32) 



-UJ 
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Parameter 



Value 
(-0.24 GeV)^ 



0.5 



(^G^) 0.012 GeV^ 
{qgG-aq)/{qq) 0.8 GeV^ 



TABLE I: Input values of the vacuum condensates at the 
normalization point p = 1 GeV. 



Cp — t^ln 
-w 

"24^^"—' 
0, 



(33) 



up to the terms polynomial in lu, which are irrelevant for 
the present purpose. Here, Cp = {N'^ — l)/(2A^c), and 
we note that the dimension-5 mixed condensate does not 
contribute to H^J'^. Substituting ([20)1 with the coeffi- 
cient functions ([32]) and dM]) into the RHS of Eqs. ^ 
and (fT9l) . respectively, we obtain 



F2(M)Al,(M)e-^/^'^: 



2a, 



^ M ' 



3a 



^^(^G^)M^(o)(^^)_l(^,G.a,),(34) 



and 



2a 



3a 



TT 

M a< 
'TV 



'^th ^ 

M ' 



(35) 



where, as usual, the factor a^/Tr from the coefficient func- 
tions Cq ' of Eqs. ((32|). ((33|) is combined with the gluon 
condensate. This set of the sum rules reproduces Grozin- 
Neubert's sum rule formulas ^ for X%{pi) and \%{p)- 

In numerical evaluations throughout this paper, we use 
the standard values for the input parameters collected in 
Table m These values have been used in, e.g., a recent 
QCD sum rule calculation for the B-meson light-cone 
distribution ainplitude [l^, and are consistent with the 
values used in Q ; note that the values of the condensates 
have been extracted at the ~ 30% level accuracy [lol |. 
as(l GeV) = 0.4 was used in 'sS], but we use the value 
as{l GeV) = 0.47 which is consistent with the world 
average. 

In Figs. [3] and m we plot X%{p-) and Xj^ifi) — A|;(^) 
a.t — 1 GeV as functions of M, obtained by tak- 
ing the ratios of Eqs. (l34)) and ([35|) . respectively, with 
Eq. (I?71) . One finds that the values of Xjj are larger than 
those of the splitting Xj^ — A|,, and, in particular, that 
the curves for the former show sizeable dependence on 
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the parameter M. Indeed, this considerable variation 
of \% for 0.3 GeV < M < 0.5 GeV, which was taken 
in 8] as the "stability window" for the sum rule (|27l) . 
is responsible for the large errors in Eq.([6|). Such poor 
stability is known to be a common feature to the sum 
rules for matrix elements of operators with high dimen- 
sion pT| . because the corresponding sum rules are dom- 
inated at small M by the condensates of high dimen- 
sion. For the present case with the dimension-5 opera- 
tors in Eqs. (Hj), ([S]), the behavior of the sum rule for 
0.3 GeV < M < 0.5 GeV is mainly determined by the 
term with the quark-gluon- mixed condensate {qgG ■ aq) : 
this is demonstrated in FigjSl which shows the separate 
contributions from each term in the RHS of (IMl) . orga- 
nized according to the dimension of the associated local 
operators. On the other hand, the term with (ggG • aq) 
is absent from the sum rule ([35]) : this sum rule yields 
a rather stable behavior for A|^ — A|; as shown in FigjH 
but the separate contributions of each term in the RHS of 
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FIG. 5: The separate contributions to the Borel sum rule 
for Xjj with = 1 GeV and ujth = l.OGeV, using Eg. pi)) 
divided by Ea. (l27p . The contributions from each term in the 
RHS of Ea. (|34|l . organized according to the dimension of the 
associated operators, are shown. 



FIG. 7: The separate contributions to the Borel sum rule 
for F^ Eq.(I13, with fi = 1 GeV, cjth = l.OGeV, and 
A — 0.55GeV. The contributions from each term in the RHS 
of Eq. (|27|l , organized according to the dimension of the asso- 
ciated operators, are shown. 
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with the dimension-5 quark-gluon-mixed condensate, is 
of 0(a°), while all the other coefficients are of 0{as). 
This is again in contrast to the case of the decay con- 
stant, for which all the nonzero coefficients in Eg. ([^5)1 
are of 0{a^). Combined with the behaviors in Figs. [5][7] 
discussed above, in particular, with those indicating the 
dominance of the term associated with the dimension-5 
operator in Ea.lfM]). it is desirable to calculate the 0{as) 
correction to this term induced by the next-to-leading 
order (NLO) correction to the Wilson coefficient C^^ 
of Eg. (15^. It should be also clarified whether C^^ of 
Eg. (155)) receive the 0{as) effects. These new 0{as) con- 
tributions may give the effects comparable with the other 
0{as) contributions displayed in Figs. [51121 We will work 
out the one-loop matching to calculate the corresponding 
Wilson coefficients at 0{as) in Sec. IIVI and derive the 
modifications to the sum rules for Xj^ ^ . 



Eg. (l35p . shown in Figl6l indicate that the nonperturba- 
tive corrections do not decrease for increasing dimension 
d = 0, 3, and 4 of the associated operators, similarly as in 
FigjS] These characteristic behaviors in Figs. [31 and HI are 
in contrast to the case of the decay constant, for which 
the separate contributions to the sum rule ([?7)) are plot- 
ted in FigjTlwith the value A = 0.55 GeV [8]. These 
results suggest good convergence of the OPE I^U^ for 
X = F, with the operators of dimension d < 5, while 
the convergence of Eg. d^ for X = 3H, 35, at the same 
level of accuracy, is questionable. Therefore, we will cal- 
culate the nonperturbative corrections to Eg. ()20p . asso- 
ciated with the dimension-6 operators, and evaluate the 
corresponding modifications to the sum rules p4l) and 
([55)) . as well as to Eq. ([?7)) . in the next section. 

The results for the relevant Wilson coefficients, 
Eqs. ([5^ and ([55)) . show that only C^^ , associated 



III. NONPERTURBATIVE CORRECTIONS 
WITH DIMENSION-6 OPERATORS 

At the leading accuracy in as, all the relevant power 
corrections to Eo. ([^n)) with X = F are generated from 
the second term in Eg. ([M)) . and the leading contribution 
in the ellipses in the middle line of Eq. ([25)) reads 

lx''x''x\0\qiO)riP+r2D^D,DxqmO), (36) 
6 

which determines the power correction associated with 
the dimension-6 operator. We exploit the following exact 
relation between matrix elements of the local operators 
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FIG. 8: A diagram for the OPE of the correlation function ([8} . 
At the leading accuracy in Os, this diagram determines the 
power correction induced by the dimension-6 operators. 



(r is arbitrary gamma matrix) [23 | , 



which can be derived straightforwardly using the equa- 
tions of motion, Ipq = 0, D^G«, = gEg'^'l^T^, 
where the summation is over all quark flavors. Thus, 
Eq.® yields 

Tr [FiP+r^] ^ [v ■ x)x' {0\qj.T^q ^ ifj^TWlO) , 

(38) 

which implies the new power-correction term in the RHS 
of Eq. (gni with X = F, given as [H 



(39) 



where, as usual, the four-quark condensate in Eq.(| 
is reduced to the square of (qq) using the factorization 
approximation through the vacuum saturation. As the 
resuh, the RHS of Eq.d^T]) receives the new term (T9|,[20t, 



72NcM^ 



(40) 



Similarly, we can calculate the new corrections to 
the sum rules ([M)) and ([55]) for A|;^, induced by the 
dimension-6 operators: at the leading accuracy in as, 
all the relevant power corrections to Eg. ([201) with X ~ 
3H, 3S are generated from the second term in Eq. p9| . 
The leading contribution in the ellipses in Eg. (1311) is as- 
sociated with the dimension-6 operators and reads 

x^{0\qiO)gG^,iO)TiP+T2DxqmO) 
iv ■ X 
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T\J2^'j''T\'\0) (Tr[a^,riP+r2 
q' 



(41) 



-|-2Tr [{iVf^^i, - Wi/7^)riP+r2]) , 



where the matrix element has been handled using Eq.([37]) 
with the indices /i and v antisymmetrized. Diagrammat- 
ically, this result is represented in FiglS) Eg. (HT]) implies 
the new power-correction term in the RHS of Eq. (PU)). 
given as 
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(42) 
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for X — 3H, 35, respectively, so that we find that 
Eqs. ([M)) . ([55]) receive the new terms. 



TrCpas /_ ,2 
2NrM ' 



TtCfQ-s i_ ,2 



(43) 



respectively, in their RHS. 

We now discuss the effect of the above-obtained power 
corrections due to the dimension-6 condensates on the 
corresponding sum rules. First of all, when including 
Eg. (HP)) in the sum rule ([TT]) . only this new term is of 
0{aa), in contrast to the other terms arising in the RHS 
of Eg. ([?7)) . and, actually, the effect of this new term turns 
out to be completely negligible. This is demonstrated in 
Fig. ini here, the dashed curve is almost indistinguishable 
from the solid curve, and the former is same as the solid 
curve in Fig. [Tj 

For the case of the sum rules ([M|) and pSI) for X\ j^, 
the terms in the RHS, except the term associated with 
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the mixed condensate {qgG ■ aq), are of 0{as) similarly 
as the new contributions of Eg.ipS]). The quantitative 
roles of these new terms are shown in Figs. [TUl [11] we 
calculate Eq. p4p with and without taking into account 
Eq.dlSl), and divide both the resuhs by Eq.dl?]), yield- 
ing the solid and dashed curves, respectively, plotted in 
Fig. [TUl The similar calculation based on Eq. (I55t yields 
the solid and dashed curves in Fig. [TT} note that the 
latter curve is same as the dot-dashed curve in Fig. [6| 
because the dimension-5 condensate does not contribute 
to Ea. (|55)) . The new contributions of Eg.lj^^]) enhance 
the values of A|j and — A|;, giving rise to some ad- 
ditional dependence on the Borel parameter M, but the 
effect is not so significant. Indeed, the comparison of 
Figs. [TU] and [TT] with Figs. [S] and [B] respectively, indi- 
cates that the effects of Eq. (|43)l due to the dimension-6 
condensates are smaller than the dominant effects from 
the lower-dimensional condensates, such that the conver- 
gence of the OPE (1201) with X ^ iS as well as X = 3i7 
may be suggested at this level of power corrections. 



ONE-LOOP WILSON COEFFICIENTS FOR 
THE DIMENSION-5 OPERATORS 



IV. 



The Wilson coefHcients arising in the OPE (pOl) are in 
general expressed as a power series in a^, 



(44) 



In particular, for the case with X = 3H, 35, the formu- 
las dSH), dag and dm indicate Cf = 0, for the 
coefficients associated with the operators of dimension 
d < 6, except for k — a and X = 3H. Those formulas 
also give the explicit nonzero results of C'l'^^^\u!) and 

3S( 

Cj, (w), except for the case with k = a, whose results 




FIG. 12: Feynman diagrams for the one-loop matching of the 
Wilson coefficients associated with the dimension-5 operators. 
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(c) 



have been unknown. In this section, we derive C, 



3H{1) 



FIG. 13: Examples of vanishing diagrams for the one-loop 
corrections to the correlator ^ in the Fock-Schwinger gauge. 



and Ca ' '(uj), performing the one-loop matching calcu- 
lation. The OPE (|20|) with this result allows us to con- 
struct the sum rules for Xjj and Xjj — X^, taking into 
account all the relevant 0{as) effects. 

The Feynman diagrams in Fig. IT^l represent the one- 
loop corrections to the correlation function in the LHS 
of Eq.®, which are relevant to the matching to derive 
Ca (i^) for X = 3H,3S; here, the corrections due to 
self-energy insertions into the quark or gluon external 
fields have been omitted, because the corresponding con- 
tributions eventually cancel in the matching. We calcu- 
late those Feynman diagrams in Z? = 4 -t- 2e dimensions 
and derive Ca (w) in the MS scheme. It is convenient 
to use the Fock-Schwinger gauge for the classical back- 
ground gluon field, as in Sec.|lTl Then, the diagrams con- 
taining the subdiagram (a) of Fig. 12] as well as the subdi- 
agram (b) or (c), vanish. Also, for the correlator the 
diagrams in Fig. 1131 vanish: in Fig. I13[ the vertex arising 
from the quark-gluon current q{0)TigG fj,u{0)hy{0) van- 
ishes using Eg. p^ . and, indeed, the corresponding con- 
tribution is absent from the relevant propagator ^U\\ . 

Here, for later use, we mention the contribution of 
the tree diagram (d) in Fig. m to the correlator in the 
LHS of Eq.®, for the case of D dimensions and under 
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the classical background fields q{x), q{x) and G^i,{x) for 
quarks and gluons; the corresponding contribution reads 
(see Eq.dni)), 



i J d^xe-'"''-''e{-v ■ x)S^^^ {x±) 

xg(0)ffG^,(0)riP+r2g(0) 

= -iTr[a^.riP+r2]nr^'=M 

-Ixr [{iv,j, - iv,j,)T,P+T2] nlTi^), (45) 

where the extraction of the relevant scalar piece from the 
arising combination of the classical fields is implicit in the 
LHS, yielding the RHS that is expressed in a similar form 
as in the RHS of Eq.®, with 



1 



-o, 



2D{D-l)uj 

nlTH = 0, (46) 

in terms of the dimension-5 quark-gluon scalar operator, 



O^qgG- aq. 



(47) 



When D 4 and O -> {qgG ■ aq), Eq.igS]) reduces to 
the corresponding contribution using the first term in the 
RHS of Eq.dSIl). 

Now, it is straightforward to calculate the diagrams 
in Fig. [T2] combining the relevant propagators ((2T|). ([23| . 
([30|) with other (familiar) building blocks, and it is conve- 
nient to perform the loop integrations in the coordinate 
space: the diagrams (e)-(g) and (i)-(j) in Fig. [T2]represent 
the contributions generated by the second term in the 
RHS of Eq.dSOl) and of Eq.®, respectively, and the dia- 
gram (h) can be calculated similarly; on the other hand, 
from the "nonlocal quark condensate" contributions con- 
tained in the eUipses in Eq. ([29l) . the diagrams (k), (1) are 
generated as the subleading terms in the Taylor expan- 
sion similar as in Eg. (1251) . We note that the diagrams (a)- 
(c) and (e) all give the UV-divergent results, while the di- 
agrams (g), (j) and (1) give the IR-divergent results. Each 
of the diagrams (f), (h), (i) and (k) vanishes as a result 
of the "canceling" UV and IR poles, l/euy — ^/^IR, aris- 
ing from the scaleless loop integral. On the other hand, 
the diagram (d) yields the result of 0(e) and vanishes as 
Z? 4. On the LHS of Eq.®, in addition to the con- 
tributions from all the diagrams in Fig. I12[ we have also 
the counter-term contribution. 



(48) 



with iqWi^) defined as Eqs. (gH), (ge]). Z,, and Z2 are 
the quark-field renormalization constants for heavy- and 
light-quarks, respectively, in the MS scheme using the 



Feynman gauge for the quantum part of the gluon field, 
as 



27re ' 



Aire ' 



(49) 



at one- loop order [2], where l/i = 1/e + jE — ln47r, with 
the Euler constant, and, similarly, Zj is the renor- 
malization constant for the heavy-light current operator 
J = hvT2q in Eq.([8|), as 



Zj 



Z., 



1 



(50) 



connecting the renormalized and bare operators at one- 
loop order 2]. In Eq.(|48|). the remaining counter terms, 
associated with the renormalization of the quark-gluon 
three-body current operator qVigG^j^^h^ in Eq.(I5]), are 
represented by the ellipses whose explicit formula is given 
in Eq. ((X3)) in Appendix El 

We combine the sum of the contributions of all the dia- 
grams in Fig. [T2] with the counter-term contribution (|48p . 
and decompose the result as in the RHS of Eq.([8]), de- 
noting the corresponding two correlation functions as 
n3"^°°P(w) and n3g°°P(a;) in place of n3H^3s(w). We ob- 
tain, for w < 0, 



n 



3H 



35 (^) 



967rw 
2 



2 

10 



- Nr - — ] In 



O, 
1 



+ 



(51) 
O, (52) 



with the ellipses denoting the terms that vanish as e — > 0. 
Here and below, fi is the MS scale and as = Q!s(/i). 
In the calculation to derive these results, we observe 
that the UV poles from the diagrams (a) and (b) are 
canceled, respectively, by the first and second counter- 
terms in Eq. pS)) . Similarly, the UV poles from the dia- 
grams (c) and (e), as well as the UV poles from the di- 
agrams (f), (h), (i) and (k) due to the above-mentioned 
l/eav — l/e/_R structure, are completely canceled by the 
UV poles arising in the ellipses in Eq.(j48|), i.e., by the 
counter terms (IA3p . Then, the remaining 1/e poles from 
the diagrams in Fig. [12] are the IR poles only, as it should 
be, and the sum of all those IR poles yields the 1/e pole 
in the above results ([?T|) . (|52p . 

As a useful cross-check of our results (ISTI) . (|5^ . we 
performed the corresponding NLO calculation also for 
the correlation function. 



dDj:e'-^--{0\T[q{x)TigGf,,{x)h,{x) h,{0)r2q{0)] |0). 

(53) 

Because the Fock-Schwinger gauge, x^Af^{x) = 0, used in 
the present paper, violates translational invariance, the 
calculation of Eq. ([5^ does not coincide on a diagram- by- 
diagram basis with that of the LHS in Eq.®; e.g., for the 
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correlator ([53]) . the contributions of the diagrams (a)-(c) 
in Fig. [T3] do not vanish. With the similar technique as 
above and with sir ~ ec/y, we calculate the contributions 
to Eq.(|53p taking into account all the relevant diagrams, 
i.e., the diagrams in Figs. [121 HSl and some other non- 
vanishing diagrams. We find that this calculation yields 
the results identical to Eqs. ([51]) . (|52|) . Indeed, for ex- 
ample, considering the diagram (e) in Fig. [T^ and the 
diagram (b) in Fig. [131 the sum of the contributions of 
these diagrams for Eq. ([53| coincides with the correspond- 
ing sum for Eq. ([8]) . This reflects the fact that the transla- 
tional invariance is restored in a gauge- invariant subset of 
the diagrams; note that the contributions of the diagrams 
involving the subdiagrams in Fig. [2] vanish for Eq. ([53| . 
as well as for Eq.0. Similarly, the diagrams (g), (j), (1) 
in Fig. [T^l the diagram (a) in Fig. [T31 and some other di- 
agrams form a gauge-invariant subset, such that they are 
obtained from the diagram (b) in Fig. [T] by attaching an 
external gluon line in all possible ways, and we observe 
that the sum of the contributions of those diagrams is 
identical between Eq.([53l) and Eq.®. With ejR = euv, 
each of the remaining diagrams in Figs. [121 [El actually 
yields the identical result for both Eqs. ([55)) and ([5]). 

The matching relations of our above results ([45)) - 
([46| . ([5T|) and ([52[) with the corresponding term in the 
OPE ^ read {X ^ 3H,3S), 



one-loop matching due to the 0{as) terms in Eq. ([54[) 
leads to the relation. 



do 



(58) 

the both sides of which are finite as e — )• 0: for X = 3H, 
the second term in the RHS serves to cancel the 1/e pole 
arising in the first term (see Eq. ([5T|) ). while, for X = 35, 

3S( D 

Eq.([57| implies that Ca ^ '(i^) is directly given by the 
coefficient of O in Eq. ([52[) . Substituting the e — ^ limit 
of Eqs. ([ST]) and ^SEIi into Eq.([ll]), we obtain the final 
form of the corresponding NLO Wilson coefficients. 
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24u} 



1 



2 Nr 



In- 



-2uj 
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967ra; 



iV.ln- 



-2w 
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2X1 



(59) 
(60) 



in the MS scheme. 



V. RENORMALIZATION-GROUP 
IMPROVEMENT AND BOREL ANALYSIS 



ntree 
X 



(^) 



n^^°°p(c.) 



(54) 



for — w ^ Aqcd, with the Wilson coefficients expressed 
as Ea. (|44p . and the renormalized composite operator 
corresponding to Eq.([47|), such that (OlO'^'lO) = 
(qgG • erg). Here, the renormalized operator 0^°^ is re- 
lated to the bare operator O^arc^ ^^^^^ |.j^g operator O 
which arises in Eqs. ([46[) . ([5T|) and ([52]) and is composed 
of the classical (renormalized) constituent fields, as 



Zo 



-O 



bare 



Zo 



(55) 



with Z2 of Eq. ([^5)) , noting that the combination gG^^i, is 
not renormalized in the background field method [23]. To 
obtain Zq , we performed the one- loop renormalization of 
the dimension-5 quark-gluon-mixed operator (j47|) in the 
present framework with the MS scheme, and the result 
is (see the discussion below Eq. (|A5[) in Appendix [5]) 



•are \ Xr 



(56) 



which coincides with the corresponding result calculated 
in Ref. Now, the matching of the 0(a°) terms of 

both sides in Eq. ([54l) immediately yields 

Cl"'^''\u) = 4(2 + g)(3 + 2e)c^' 

Cf(°)(a.)=0, (57) 

for arbitrary e, and these results reproduce the corre- 
sponding formulas in Eqs. I\'62\i and ([55)) as e ^ 0. The 



We substitute the new results ([59)) and 
Wilson coefficient (oj) of Eq.([20l), while. 



([60)) into the 
for the other 



coefficients Cf (w), (w), and Cq (w), we have the cor- 
responding formulas in Eqs. ([5^ and ([55)) : furthermore, 
we add the new power-correction term with (|42p to the 
RHS of Eq.([20l). The result gives our upgraded OPEs 
for the correlation functions in Eq.([51), taking into ac- 
count the operators of dimension d < 6 and the asso- 
ciated Wilson coefficients to the 0{as) accuracy. Now, 
we use these new results for the OPE to derive the sum 
rules for Xj^ and Xj^ — X^' we substitute these OPEs 
into Eqs. ([TSIl and ([TO)) : here, it is straightforward to 
calculate the discontinuities of the coefficients ([59)) and 
([60| across the cut along the line w > in the complex 
ll) plane, reexpressing the logarithmic contributions as 
{2 / u})\n{—2u} / fj,) = {d / du))\v? {—2u) / ^) . As a result, we 
obtain the new formulas of the Borel sum rules for A|j 
and A|j — A|; , 
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-{qgG-aq}-^ 
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2M , Wth ^ 
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2^. 



nCpOls 



where T^("'(wthA/) is defined as Eq.dSH]) and 



r(a, z) 



(62) 



(63) 



is the incomplete Gamma function. These two sum 
rules, (pT|) and (15^ . are the new results that take into 
account the operators of dimension d < 6 and the 
associated Wilson coefficients to the 0{as) accuracy: 
compared with the previous results ([M)) and (1551) that 
correspond to Grozin-Neubert's sum rule formulas 8], 
Eqs. ([6T|) and (|62l) receive the 0{as) corrections associ- 
ated with the dimension-5 quark-gluon-mixed condensate 
(qgG-aq), as well as Ea. (|i5)) due to the dimension-6 four- 
quark condensate (99)^ • In particular, the former correc- 
tions bring an explicit dependence on the scale /i to the 
RHS of Eqs. (|5T|) and ([5^. through the logarithmic term, 
ln(2M/ fie'"^): one can show that (m) determined by 
our formulas (j6ip and (j62p satisfy the renormalization- 
group equations of Eqs. (jA4l) . (lASp . taking into account 
the derivative of the above logarithm \n{2M/ ne""^), as 
well as the scale dependence of the other terms controlled 
by the nontrivial anomalous dimensions: 



(64) 
(65) 
(66) 



M^+7./(as) )F(^)=0, 



+7g(as) ) {qq)if^) = 0, 



7fc(as) = 7fcO^ +7fci (^j H (fc = J, q, cr), (67) 

where 



7J0 = -3Cf, 790 = -6Cf, 7cto ^ N^- ^,(68) 



7ji = Cf 



5 _ 8^ 
2 3' 
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7,1 = -Cf ( 3Cf + y - —Ni 



(69) 
(70) 



while 7cti is not available. For the present purpose to 
confirm Eq. (|A4p with the one- loop mixing matrix (|A5p . 
the explicit forms of 7jo and 7cto in Eg. (|68l) . combined 
with dk/dn = 0, are sufficient, where the former 
are immediate consequences of the corresponding renor- 
malization constants ([5(I| . (|56l) . Apparently, Grozin- 
Neubert's sum-rule formulas given by Eqs. and ([55)) 
do not obey such renormalization-group property. Thus, 
our sum rules (1611) and (I62p allow the first nonpertur- 
bative estimate of the HQET parameters \\ with 
the correct /i-dependence implemented. We emphasize 
that the new 0{as) corrections in Eqs. (IFTI) and (|62p . 
associated with the dimension-5 quark-gluon-mixed con- 
densate l^gG ■ aq), play essential roles to reproduce 
renormalization-group equations of Eqs. (IA4p . (IA5p in 
the QGD sum-rule framework. 

It is worth comparing this remarkable property of 
Eqs. (pT|) and ([5^ with the situation for the case of the 
sum rules of the decay constant F{^), based on the cor- 
relator ([7]): F{fi) determined from the sum rule ([27]), 
which was presented in Sec. [iTj at the leading accuracy 
in as, does not obey Eq. ([M|) . even though we take 
into account the scale dependence of Eqs. ([S5|) . (|66l) in 
Ea. (l27p . Now, including the higher-order contributions, 
the 0{as) corrections for the relevant Wilson coeffi- 
cients and the dimension-6 condensate contribution (j40p , 
Eq.([27| is modified into i, |M HI 
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Here, in particular, the 0{as) corrections arising in the 
second line bring an explicit dependence on the scale /i, 
through the logarithm ln(/i/2Mz). Taking into account 
this new /i-dependence, F{fi) determined by Eq. ()7ip 
obeys the renormalization-group equation (|64p , up to the 
corrections of 0{aj.) and the small contributions from the 
condensates of dimension-5 and higher (see Figs. [3 (HI). 

The correct renormalization-group properties obeyed 
by Eqs. ([CT|) . ([5^ and ([7T|) allow us to improve these 
sum rules further, such that the logarithmic effects asso- 
ciated with Q;sln(M/ //) are resummed to all orders. This 
renormalization-group improvement is formally achieved 
by setting /i = /i' - M in Eqs. ([H]), ([H]) and ([H]), fol- 
lowed by evolving the resulting HQET parameters and 
the condensates at the scale ^j.' to those at /x ~ 1 GeV: 
using the first two coefficients of the (3 function, 
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with Nf being the number of active flavors, the corre- 
sponding renormahzation-group-improved sum rule for 
the decay constant reads 
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Here, we have neglected the unknown NLO-level effects 
associated with the dimension-5 quark-gluon-mixed con- 
densate, i.e., the corresponding one- loop coefficient func- 
tion and two-loop anomalous dimension, and also ne- 
glected the running of the dimension-6 four-quark con- 
densate. Up to these small effects, Eq.(|73| combined 
with Ea. (l74l) sums up the leading and next-to-leading 
logarithms of the form <ln" (M//i) and <+iln"(M/^). 

As a result of the similar renormalization-group im- 
provement for the two sum rules (j6ip and (j62p . we obtain 
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Here, the anomalous dimensions for //(m) are known 
only at one- loop as Eq. ljASp . We can divide Eqs. (|75|) 
and ([75)1 by Eq.([75| so as to eliminate the factors 
F2(^')e-'^/*^ before comb ining with the evolution ((77|). 

We evaluate the renormalization-group-improved 
Borel sum rules Eqs. ([731), ^ and ([751) with fi' = 2M 
and fj, = 1 GeV; here, the choice /i' = 2M allows us 
to sum up the relevant logarithmic contributions. In 
addition to this default choice, we also study the scale 
dependence of the results when varying fi' in a range 
with /i' ~ M; this gives us a rough estimate of the 
uncertainty of the results due to the neglected part of 
order a^, and to the higher-order contributions. We use 
the input values for the vacuum condensates as given in 
Table |T] in Sees. [TTl and use the two-loop expression for 
the running coupling as(/i) with Aq^q = 0.31 GeV, so 
that as{l GeV) ~ 0.47, and as(mB) ~ 0.21. 

First of all, we note that 
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corresponding to the /j,'-dependent factors in Eq.lfM]). 
forms the renormalization-group-invariant combination, 
for which the Borel analysis has been performed in the 
literature 0,[ll|23: substituting Eq.dTS]) into the RHS 
of Eq. l|75|) . we obtain the sum rule formula for the 
renormalization-group-invariant decay constant. 
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FIG. 14: Borel sum rule for A based on Eg. (TfO)) with ^J,' 
2M. From bottom to top, the curves correspond to ujth 
1.1,1.2,1.3,1.4,1.5,1.6 GeV. 



FIG. 15: Borel sum rule for F based on Eq.(l79| with ^' = 2M 
and A = 0.4 GeV. From bottom to top, the curves correspond 
to tJth = 1.1, 1.2, 1.3, 1.4, 1.5, 1.6 GeV. 
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where ^' = 2M as the default scale. Then, taking the 
ratio of Ea. ((75|) and its first derivative with respect to 
1/M yields an equation to estimate the value of A. Us- 
ing this equation, Fig. [14] displays A as a function of the 
Borel parameter M for different values of the threshold 
parameter wth- The stability of this sum rule is poor 
for M < 0.5 GeV because of the growth of as{2M), 
and such behavior is more pronounced in the present 
case with as{l GeV) ~ 0.47 than in the literature M 
with as(l GeV) ~ 0.34 (see also the discussion in [2|. [20j 
to avoid the popular choice fi' = 2M). Also, our re- 
sults suggest smaller values for A and larger values for 
wth than those in 0, H, As a result. Fig. [T3] shows 
that the values 0.4 GeV < A < 0.5 GeV correspond to 
1.2 GeV < wth ^ 1-5 GeV. Indeed, the perturbative cor- 
rections tend to increase the pole mass [l^, [2^, [27j 
and thus tend to decrease A = tub — mi,; however, 
A < 0.3 GeV corresponding to Wth < 1-1 GeV in Fig. [Til 
would be somewhat too small in view of the experimental 
information on the 6-quark mass (28j . 

In Figs. [T5| and [TBI we show F as a, function of M 
using Eq.((7n|) with A = 0.4 and 0.5 GeV, respectively, 
where the curves are drawn for different values of Wth- 
The stable behaviors are obtained for M > 0.4 GeV with 
1.2 GeV < wth < 1.4 GeV, yielding F ~ 0.4 GeV^/^. 
It is worth noting that, by taking into account the cor- 
rections due to the finite quark mass nib, the value of 
F around 0.4 GeV^/^ was shown to be modified into 
Jb — 0.2 GeV 0, H, which appears to be consis- 
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tent with the recent lattice results, Eq.Q. We mention 
that Eq.([79l) with A > 0.5 GeV yields the larger F as the 
stable value associated with wth ^1-4 GeV; this would 
not be favored in view of Eq.Q. Thus, we conclude that 
the range wth — 1.2—1.4 GcV corresponds to the opti- 
mal choice for the threshold parameter, and, with this 
choice, we will evaluate Borel sum rules for ^; for 
this purpose, we need not specify the value of A because 
Eqs. ([75)1 and (1751) divided by Eg. ([75)1 are independent 
of A. By inspection of the contributions from each term 
in Eq.dfSl), we find that, for M > 0.4 GeV, the con- 
tribution of the perturbative correction terms are less 
than ^ 40% in the OPE for the sum rule and the con- 
tribution of the nonperturbative power correction terms 
is much smaller. On the other hand, for M < 0.6 GeV, 
the contribution of the higher resonances and continuum 
contributions is less than ~ 40% of the total contribu- 
tion in the dispersion relation for the sum rule. Thus, we 
take 0.4 GeV < M < 0.6 GeV as the stability window in 
the following calculations of Eqs. ([75]) and ([75)) . Indeed, 
one can also confirm that the correction terms, as well 
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FIG. 18: Same as Fig. [13 but for A|(1 GeV). 



as the higher resonances and continuum contributions, 
arising in the sum rules (|75p . (|76p. are under control for 
0.4 GeV < M < 0.6 GeV. 

In Figs. [T7] and [TH] we shovif the results for A|^(/i ~ 
1 GeV) and — 1 GeV) , respectively, as functions of 

the Borel parameter M , using the continuum threshold 
Wth = 1.3 GeV. The solid curves are obtained by dividing 
Eqs. ([7S|) and ^ by Eq.([731), followed by the evolution 
of the results from /x' = 2M to ^ = 1 GeV using Eg.dTTl). 
and thus show our full results with the renormalization- 
group improvement. For comparison, the dot-dashed 
curves show the "fixed-order" results that are obtained by 
dividing Eqs. ^ and ([62]) by Eq.([71]), while the dashed 
curves are obtained by dividing Eqs. (j34[) and psp by 
Eq. (|27| . The dashed curve in Fig. [17] shows the behav- 
ior similar as the curves in Fig. [3] because the differ- 
ence between those curves are only due to the value of 
the threshold wth. The new 0{as) contributions due to 
Fig. \12\ and also the associated renormalization-group- 
improvement effects, significantly improve the stability 
of the sum rules. Furthermore, those contributions sig- 
nificantly reduce the magnitude of as well as \\. For 
those results, we find that the behavior of the decay con- 
stant arising in Eqs. (|75)) and ([751) . actually plays im- 



portant roles: comparing Fig. [9] with Figs. [15] [16] we 
see the large effects due to the renormalization-group- 
improved, NLO perturbative corrections in the decay 
constant sum rule. It is worth mentioning that the corre- 
sponding large radiative corrections to the decay constant 
is mainly due to one-gluon exchange between heavy and 
light quarks in Feynman gauge, i.e., from their Coulomb 
interaction, and that those large effects are essential to 
yield the values consistent with Eq.® 0,0; on the other 
hand, it has been demonstrated that the decay constant 
sum rule (|73p is quite stable with respect to inclusion of 
the NNLO-level radiative corrections [lOj . 

Figs. [T9l and [20l show the ujth dependence of the Borel 
sum rules for A|^(l GeV) and A|;(l GeV), respectively, as 
functions of A/, which are obtained by dividing Eqs. ([75)1 
and ([7S|) by Eq. ([75)) . followed by the evolution of the 
resuhs from ^' = 2M to ^ ^ 1 GeV using Eq. ([77)) : 
the middle curves in Figs. [TO] and [201 are same as the 
solid curves in Figs. [T7] and [TH] respectively. We see 
that the central values are A|j(l GeV) = 0.06 GeV^ and 
A|;(l GeV) — 0.03 GeV^ as an average over the sta- 
bility window M = 0.4—0.6 GeV and the range wth = 
1.2—1.4 GeV for our optimized threshold parameter, and 
that these central values are associated with the uncer- 
tainties ±0.01 GeV^ and ±0.015 GeV^, respectively. 

To estimate the uncertainties due to the lack of in- 
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formation of two-loop anomalous dimensions as well as 
the higher- loop effects in Eqs. ([75|) - d771) . the scale depen- 
dence of our full results with the renormalization-group 
improvement is analyzed by varying around the de- 
fault value 2M, e.g., for a range M < ^' < AM; if 
fi' — M were chosen, as{M) would arise in the cor- 
responding formulas with too small scale for perturba- 
tion theory in 0.4 GeV < M < 0.6 GeV. Thus, we 
vary ^' in the range 1.5 A/ < ^' < 4M, for Eqs. (|75)) 
and ([751) divided by Ea.([75|. followed by the evolution 
using Ea. ([77|) from ^' to ^ = 1 GeV: the correspond- 
ing results of A|; ^(1 GeV) increase (decrease) for in- 
creasing (decreasing) /i', and we find that the above- 
mentioned central values, A|^(l GeV) = 0.06 GeV^ and 
AKl GcV) = 0.03 GeV^, receive the ±0.02 GeV^ and 
±0.005 GeV'^ variations, respectively. 

Similarly, we calculate the uncertainties of the results 
due to the uncertainties in the input parameters of Ta- 
ble HI which are known to be ~ 30% or less 10]. We also 
vary A^^^^ in a range 0.29 GeV < A^^^^ < 0.33 GeV cor- 
responding to as{l GeV) = 0.44—0.5. Among them, the 
uncertainty of the dimension-5 quark-gluon-mixed con- 
densate {qgG ■ crq) produces the largest effect as ~ 15% 
and ~ 30% of the total contribution to A|^(l GeV) and 
A|;(l GeV), respectively, while each of the other uncer- 
tainties yields 10% or less of the total contribution to 
A|^^(l GeV). 

Adding the errors induced by all source of uncertain- 
ties, wth, ^J■\ condensates in Table HI and Aq^j-,, dis- 
cussed so far in quadrature, we obtain ±0.025 GeV^ and 
±0.018 GeV^ for X%{1 GeV) and A|(l GeV), respec- 
tively. There also exists an overall intrinsic uncertainty 
of the QCD sum rule method itself which is difhcult to 
estimate. Thus, with a conservative estimate of the un- 
certainties, our final results read 

A|(l GeV) = 0.03 ± 0.02 GeV^ 

A|(l GeV) 0.06 ± 0.03 GeV^ (80) 

Note that the errors in the previous estimate (|6]) Q are 
due only to the choice of the continuum threshold and the 
dependence on the Borel parameter in the corresponding 
sum rules using Eqs. (p4| and ([35]) . 

VI. CONCLUSIONS 

We have discussed the QCD sum rule calculation of the 
HQET parameters and Xjj, which represent quark- 
gluon three-body components in the _B-meson wavefunc- 
tion. We have updated the sum rules for ^ calcu- 
lating the new higher-order contributions to the OPE 
for the corresponding correlator, i.e., the order as radia- 
tive corrections to the Wilson coefficients associated with 
the dimension-5 quark-gluon mixed condensate, and the 
power corrections due to the dimension-6 vacuum con- 
densates. Combining with the similar NLO-level calcula- 
tion for the decay-constant sum rule which is consistent 



with the precise result from recent lattice QCD calcula- 
tions, we have constructed the Borel sum rules for ^. 
We have found that the new order-as radiative correc- 
tions significantly reduce the values of X% ^1^° 
make the corresponding sum rule formulas for ^ obey 
the correct renormalization-group equations. The resum- 
mation of the relevant logarithms of the 6-quark mass 
based on the renormalization group has been performed 
and proves to improving the stability of the correspond- 
ing Borel sum rules. Our final results are obtained as 
Eq. ()80|) . where the perturbative as well as nonperturba- 
tive corrections are under control and the various sources 
of errors are taken into account. 

Compared with the previous estimate, Eq.®, ob- 
tained in the central values of our results ([80)) are 
smaller by 1/3 and the errors are also reduced consider- 
ably. On the other hand, the upper bounds of our re- 
sults (|5D|) are close to the lower bounds of Eq. © . Study 
of the _B-meson light-cone distribution amplitudes using 
the new results in the present paper will be presented 
elsewhere (33l |. 
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Appendix A: one-loop renormalization of 
dimension-5 heavy-light operators 

In this Appendix, we calculate the one-loop corrections 
for the quark-gluon current operator arising in Eq.(|8l), 

qgG^iyTihv, (Al) 

in the HQET in _D 4 + 2e dimensions, and determine 
the corresponding renormalization constants in the MS 
scheme. Based on this result, we also write down the 
explicit formula implied by the ellipses in the counter- 
term contribution (|48p . 

In Eq. ljAip . Fi is an arbitrary gamma matrix, and 
/i, v are the free Lorentz indices. The relevant one- 
loop diagrams are shown in Fig. 1211 and those loop cor- 
rections for the operator (jAip induce the mixing with 
the dimension-5 heavy-light operators having the same 
Lorentz-transformation property as that of Eq. (|Aip . The 
corresponding mixing matrix can be obtained, in princi- 
ple, by calculating only the one-particle irreducible dia- 
grams (a)-(h) in Fig. [5T1 but it is known, for the case 
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FIG. 21: The Feynman diagrams relevant for the one-loop 
renormalization of the quark-gluon three-body operator (lAlfl . 



of the higher dimensional operators like Eg. dAll) . that 
such procedure requires to treat explicitly the additional 
mixing of the operators that would vanish by the use of 
the equations of motion [29j : in particular, this implies 
that we have to use the off-shell external fields which, 
in turn, allow the mixing of the gauge-noninvariant op- 
erators as well as the gauge-invariant ones, via the cor- 
responding one-particlejrreducible diagrams, as demon- 
strated in many works 30] for the case of the renormaliza- 
tion of the higher twist operators relevant to the nucleon 
structure functions. 

Here, to avoid such complication associated with the 
mixing of those "alien" operators, we calculate the rel- 
evant one- particle-reducible diagrams (i), (j) in Fig. [2T] 
as well, so that we can adopt the usual background field 
method, as in the calculations of the correlators discussed 
in the main text; then, we may use the equations of mo- 
tion for the external fields at any step of calculation (3ll | 
and the contribution from each diagram of Fig. [5T]is ob- 
tained in a gauge-invariant form expressed solely in terms 
of the (many) three-body operators of the similar type 
as Eq.(jAl|. We use the building blocks ([2T1) . (|23l) . (pCT| . 
etc., and carry out the loop integrations in the coordi- 
nate space. The contributions of the diagrams (f) and 
(j) in Fig. \n] vanish because these diagrams contain the 
vanishing subdiagrams (a) and (b) in Fig. [51 respectively. 
Also, the diagrams (a) and (b) in Fig. [5T] vanish by the 
reason similar as the diagrams in Fig. 1131 The sum of the 
UV poles from all the other diagrams in Fig. [21] reads, 
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with the last line indicating that all the preceding terms 
have to be antisymmetrized under the interchange fi ^ v. 
Here, the first three lines correspond to the contributions 
of the diagrams (c), (d), and (e) in Fig. [211 the fourth 
line corresponds to the diagrams (g) and (h), and the 
fifth and sixth lines correspond to the diagram (i). 

Now, the counter-term contribution for the renormal- 
ization of the operator (|A1[) should be combined with 
Ea. (IA2l) . so as to cancel all the UV poles arising in 
Ea. (IA2l) . Thus, the corresponding counter-term contri- 
bution is given by the minus of Ea. (|A2p with the re- 
placement, 1/e 1/e (= 1/e -|- 7_b - ln47r), in the MS 
scheme. Substituting this counter-term contribution for 
the quark-gluon current (IA1[) into Eq.® and evaluat- 
ing the corresponding correlator in D dimensions under 
the background fields, we immediately find that the for- 
mula implied by the ellipses in Eq. (j48|) is formally given 
by Eg. (1451) with q{0)gGpi,{0) replaced by the minus of 
Ea. (IA2l) with hy ^ I, e ^ i, and reads, using Ill^^''{uj) 
of Eq.®, 
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where the final form in the RHS is presented up to the 
irrelevant terms that vanish as e — 0. 

It is worth mentioning that the results based on 



Eq. (|A2|) with the particular choices, Fi = and 
Fi — 7^75?;'^, reproduce the one-loop renormalization 
mixing matrix between the two pseudoscalar operators, 
qgG ■ a^^hy and qgGf^yj^v'^jc^hy, calculated in In 
particular, it is straightforward to see that the matrix el- 
ements of a system of the corresponding operator-mixing 
formulas yield the one-loop renormalization group equa- 
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tions for A| ^ of Eqs. dH, (P, as [32 



d / A|(/i) \ Osi^J^ ( A|(/i) 
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with the mixing matrix, 

\Cf + In, 



70 = I 



= 0, 



(A4) 



(A5) 



Finally, we mention that the one-loop renormalization 
of the operator (|47p to obtain the renormalization con- 



stant (|56p can be carried out in a similar manner as 
above. Indeed, most results can be obtained from the 
above result (jA2|) by the formal substitutions, hy q 
and Fi — >■ cr^^, for an external the quark field and the 
gamma matrix structure, respectively, and thus need not 
new calculation. The only diagram that requires new 
calculation is the one corresponding to the diagram (e) 
in Fig. [211 We have only one flavor-singlet, scalar three- 
body operator of dimension 5, given by Ea. (|T7)) . and thus 
do not encounter the operator mixing in the background 
field method for this case. 
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